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ABSTRACT 

^  In  actual  practice ,  iteration  methods  applied  to  the  solution  of  finite  systems  of 
equations  yield  inconclusive  results  as  to  the  existence  or  nonexistence  of  solutions 
and  the  accuracy  of  any  approximate  solutions  obtained.  On  the  other  hand,  construction 
of  interval  extensions  of  ordinary  iteration  operators  permits  one  to  carry  out  interval 
iteration  computationally,  with  results  which  can  give  rigorous  guarantees  of  existence 
or  nonexistence  of  solutions,  and  error  bounds  for  approximate  solutions.  Examples  are 
given  of  the  solution  of  a  nonlinear  system  of  equations  and  the  calculation  of  eigen¬ 
values  and  eigenvectors  of  a  matrix  by  interval  iteration.  Several  ways  to  obtain  lower 
and  upper  bounds  for  eigenvalues  are  given.  J- 
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SIGNIFICANCE  AND  EXPLANATION 


order  to  complete  the  solution  of  a  problem  in  applied  mathematics  and  obtain 
t:i  ■  results  useful  form,  one  must  often  solve  a  system  of  linear  or  nonlinear  equa- 
ti  ,  or  compute  eigenvalues  and  eigenvectors  of  a  matrix,  which  is  a  special  case  of 
.1  onl inear  system  of  equations.  A  frequently  used  method  for  solution  of  equations 
iteration,  which  has  been  investigated  extensively  from  a  theoretical  standpoint, 

■u.d  jften  applied  successfully.  In  order  to  use  iteration,  however,  one  must  have  an 
initial  approximation  to  the  solution;  furthermore,  even  if  the  results  obtained  look 
good ,  the  questions  of  existence  of  a  solution  and  accuracy  of  the  numerical  solution 
are  not  conclusively  answered.  By  use  of  interval  iteration,  however,  employing  what 
are  called  computable  interval  extensions  of  ordinary  iteration  operators,  it  is  pos¬ 
sible  to  examine  regions  for  existence  or  nonexistence  of  solutions,  and  carry  out 
the  iteration  in  such  a  way  that  lower  and  upper  bounds  for  the  solutions  are  given 
at  each  step  and  in  the  final  result  on  the  basis  of  the  actual  computation.  Thus, 
questions  of  existence  and  accuracy  can  be  answered  without  tedious  analysis,  and  one 
can  have  confidence  in  the  answers  obtained  by  interval  computation.  Now  that  it  is 
possible  to  perform  interval  arithmetic,  for  example,  at  speeds  comparable  to  ordinary 
arithmetic,  interval  solution  of  systems  of  equations  becomes  ible  method  in  prob¬ 

lems  in  which  accuracy  is  important,  such  as  the  design  of  struetur.  airframes,  and 


summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


SOLUTION  OP  FINITE  SYSTEMS  OF  EQUATIONS  BY  INTERVAL  ITERATION 


L.  B.  Rail 


1 .  Systems  of  equations.  Systems  of  v  equations  in  \>  unknowns  are  often  presented 
in  the  form 

(1.1)  f (x)  -  0 

V  V 

for  a  zuio  x  -  x*  of  the  operator  f:  D  C  R  -+  R  ,  or  as 

(1.2)  *  -  $(x) 

for  a  fixed  point  x  ■>  y  *  of  $ :  0  C  ru  -♦r'j.  The  latter  formulation  will  be  considered 
herei  there  are  many  ways  to  transform  (1.1)  into  an  equivalent  fixed  point  problem  (1.2). 
An  example  of  a  transformation  of  this  type  will  be  given  below. 

The  form  of  (1.2)  suggests  the  frequently-used  method  of  iteAAtiont  Starting  from 
a  point  xQ ,  in  the  belief  that  it  is  a  more  or  less  good  approximation  to  x*,  the  se¬ 
quence  (xn)  is  generated  by 

(1.3)  xn+1  -  $<x  ),  n  “  0,1,2,...  . 


In  case  (1.1)  has  been  transformed  into  (1.2),  $  is  often  called  an  itZAation  opeAato/i 
for  the  solution  of  (1.1).  In  the  following,  it  will  be  assumed  that  f,$  and  other  op¬ 
erators  considered  are  continuous.  Under  this  assumption,  it  follows  that  the  convergence 
of  the  sequence  {xr)  implies  the  existence  of  a  fixed  point  x*  of  if,  that  is, 


lim  ,  . 

(x  ) 
n-«°  n 


satisfies  (1.2).  (Of  course,  if  if  does  not  have  a  fixed  point,  then  {x^ )  cannot  converge; 
thus,  nonexistence  of  a  fixed  point  x*  of  if  implies  divergence  of  {*„}•) 

This  is  all  very  well,  but  in  actual  practice,  one  can  neither  represent  real  num¬ 
bers  nor  perform  real  transformations  exactly;  instead  of  {xn),  one  gets  a  finite  sequence 
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N 

(z  }  „  such  that  z„  •  x„ ,  but  z  F  x  otherwise.  In  order  for  {z  )  to  be  useful, 

n  n-o  o  u  n  n  n 

one  ordinarily  has  to  show  that  {x  }  converges,  and  then  produce  comparisons  of  x  with 

n  n 

x*  and  z  with  x  ,  a  problem  which  can  be  difficult  in  itself,  or  at  least  tedious.  A 
n  n' 

more  straightforward  computational  approach  to  the  approximate  solution  of  (1.2),  based 
on  the  methods  of  interval  analysis,  will  be  given  in  what  follows. 

2.  Interval  analysis.  Just  as  real  analysis  is  concerned  with  the  study  of  trans¬ 
formations  of  real  numbers  and  vectors,  interval  analysis  [9],  [11]  deals  with  the 
same  for  nonempty  finite  closed  intervals 

(2.1)  X  »  [a,b]  »  (x  |  a  s  x  s  b,  x  £  R) 

of  real  numbers,  and  n-tuples  of  such  intervals,  called  inteival  vectOXi.  Interval 
analysis  is  related  to  real  analysis  in  two  important  ways.  First,  the  real  numbers 
can  be  identified  with  a  subset  of  the  set  of  real  intervals  (2.1),  namely,  the  set  of 
degenerate  intervals  with  equal  endpoints.  For  a  real  number  x,  this  identification  is 
symbolized  by 

(2.2)  x  -  (x,x) . 

Second,  real  transformations  have  exte.ni.ioni  to  interval  transformations,  as  will 
be  discussed  in  more  detail  below,  and  thus  can  be  considered  to  be  restrictions  of  in¬ 
terval  transformations.  In  spite  of  this,  interval  analysis  does  not  subsume  real  anal¬ 
ysis,  but,  like  complex  analysis,  is  a  distinct  branch  of  mathematics  with  its  own 
theory,  techniques,  and  applications. 

Strictly  speaking,  an  interoat  exteniion  t  of  a  real  transformation  ♦  has  the  prop¬ 
erties  of  inclusion 

(2.3)  *(x)  -  (♦(*)  1  x  e  x)  c  $(x) 
and  restriction 

(2.4)  4(x)  -  <Mx) 

in  the  sense  of  (2.2).  An  important  property  of  certain  interval  transformations  is 
inclusion  monotonicity:  $  is  said  to  be  inclusion  monotone  if 
xcpt(x)  C*(2), 

-  2  - 


(2.5) 


Interval  extensions  can  be  formed  of  the  operations  of  real  arithmetic,  considered 


2 

as  transformations  from  R  into  R,  that  is,  of  f(x,y)  *  x  +  y,  etc.,  with  the  results 
forming  the  uxtVLoal  a/iithmztic  of  R.  E.  Moore  [9],  [11].  For  example,  muttiptication 
of  intervals  is  defined  by 

(2.6)  [a,b]*[c,d]  =  (min  IT,  max  IT],  n  *  (ac ,ad ,bc ,bd} , 

and  so  on.  These  extei  sions  are  inclusion  monotone,  and  can  be  used  to  form  inclusion 
monotone  interval  extensions  of  rational  real  functions  automatically,  simply  by  re¬ 
placing  the  variables  involved  by  interval*,  and  the  arithmetic  operations  by  their 
interval  counterparts  [9],  [11]. 

A  special  class  of  intervals  which  will  be  useful  later  is  the  set  of  'j yrmzXAAjC. 
intervals,  which  are  intervals  S  of  the  form 

(2.7)  S  =*  [-s,sj  -  s-  [-1,1] ,  s  >  0. 

More  generally,  interval  vectors  and  matrices  are  said  to  be  symmetric  if  their  com¬ 
ponents  are  intervals  of  the  form  (2.7).  It  follows  from  the  definition  of  interval 
addition  I  9J/  fllj, 

(2.8)  [a,b]  +  [c,d]  »  (a  +  c,b  +  d] 

that  the  sum  of  symmetric  intervals  is  symmetric.  Furthermore,  let 

(2.9)  | [a,b] |  =  max{ |a| , |b| }. 

Then ,  from  (2.6)  , 

(2.10)  [a,b]  •  t-s,s]  - 

hence,  multiplication  of  an  arbitrary  interval  by  a  symmetric  interval  gives  a  sym¬ 
metric  interval. 

Furthermore,  for  x  =  [a,b] ,  m(X)  *  m((a,b])  *  (a  +  b)/2,  w(X)  =  w([a,b])  *  b  -  a, 
one  cai»  write 

(2.11)  X  =  m  (X)  +  j*r(X)  •  [-1,11 

in  terms  of  the  midpoint  m(X)  and  width  w(X)  of  the  interval  X.  The  representation 

(2.11)  of  the  interval  x  is  sometimes  called  its  midpoint -hat (width  or  centzAed  form. 
Multiplication  of  X  by  a  real  number  a  can  thus  be  expressed  by 
«*X  =  tin  (X)  +  ji  i|'w(X)  •  [-1,1]  • 
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(2.12) 


The  definitions  of  |x|,  m(x),  and  w(x)  extend  to  interval  vectors  and  matrices 
by  componentwise  interpretation,  for  example,  m(X)  “  (m(X^)  , . . .  ,m(X  ) )  ,  and  so  on. 

Formula  (2.11)  then  holds  without  modification.  Formula  (2.12)  extends  immediately 

to  the  inntA  product  (a,x)  of  a  real  vector  a  •  (a1,...,a  )  with  an  interval  vector 

X  *  (x, , . . . ,x  ) .  One  has 
1  '  v 

V  , 

(2.13)  <a,X)  »  l  o-X.  »  (a,m(X) )  +  -( | a| ,w (x) ) •  [-1 ,1]  , 

i-1  1  1 

where ,  of  course ,  | a |  ■  ( | | , _ _ | | ) . 

Interval  extensions  can  also  be  made  of  other  than  rational  functions,  and  such 
extensions  can  be  combined  with  interval  arithmetic  to  extend  a  wide  claBS  of  real 
functions.  A  useful  technique  for  obtaining  interval  extensions  of  smooth  real  func¬ 
tions  is  based  on  the  mean  vaZue  (0*Jr>:  Suppose  that  <j>  is  (Frdchet)  differentiable, 
and  its  derivative  $ 1  has  the  inclusion  monotone  interval  extension  $ 1 .  Then,  t  de¬ 
fined  by 

(2.14)  *(X)  -  *(y)  +  «'(X)*(X  -  y),  y  -  m(X)  , 

is  an  inclusion  monotone  interval  extension  of  $  (  3 ] »  where  the  indicated  matrix- 
vector  operations  are  performed  in  interval  arithmetic.  Using  (2.11),  (2.14)  can 
be  written 

(2.15)  *(X)  -  ♦(y)  ♦!«•  «>•««)•  [-1.1J,  y-m(x), 

so  the  mean  value  form  t (X)  can  be  expressed  as  the  sum  of  the  real  vector  ♦ (y)  and 
a  symmetric  interval  vector.  This  elementary  observation  will  be  useful  later. 

3.  Interval  computation .  In  actual  practice,  computing  is  done  using  a  finite 
set  G  of  rational  numbers,  rather  than  the  set  of  real  numbers  R.  Interval  analysis 
can  be  adapted  readily  to  G,  which  is  also  an  ordered  set.  Extension  of  the  results 
to  be  described  here  for  one  dimension  can  be  made  in  a  componentwise  fashion  to  in¬ 
terval  vectors  and  matrices.  The  key  to  interval  computation  using  only  elements  of 
G  is  the  operation  of  diAecfed  tcundirrg,  which  will  now  be  explained. 
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Let 


(3.1)  2.  3  |  g  €  g)  ,  g  *  max{g  |  g  £  G} . 

Attention  will  be  confined  to  the  set  of  real  numbers 

(3.2)  RG  *  {x  |  i  S  x  S  g,  x  6  R}, 

and  attempts  to  compute  a  real  number  z  such  that  z  <  j  or  z  >  g  will  be  said  to  re¬ 
sult  in  an  error  condition  called  ov&ifiZow.  (Alternatively,  the  extended  real  numbers 
could  be  adjoined  to  G,  with  overflow  detected  in  case  the  result  of  a  computation 
is  a  semi-infinite  or  infinite  interval.)  Now,  IG  will  denote  the  set  of  all  inter¬ 
vals  with  endpoints  in  G,  that  is, 

(3.3)  IG  -  { [a,b]  |  a  s  b,  a,b  £  G}. 

This  is  the  set  of  intervals  which  are  txaztly  HipHej,&niabt&  in  terms  of  the  numbers 
which  are  available  in  the  actual  computation.  For  x  £  rg,  the  upwa/id  Hounding  op&H- 
aXoH  A  to  G  is  defined  by 

(3.4)  Ax  ■  min{g  |  g  2  x,  g  £  G) , 

and  the  doujnuxvid  Hounding  openatOH  7  to  G  by  the  corresponding  expression, 

(3.5)  Vx  ■  maxfg  |  g  s  X,  g  E  G) 

(  7  1  .  For  the  set  IRG  of  intervals  with  endpoints  in  RG,  application  of  the  diHictid 
Hounding  optAaiOH  □  to  intervals  [a,b]  £  irg  gives 

(3.6)  D[a,b]  -  [7a, Ab]  £  IG, 

and  thus  □  maps  IRG  into  IG.  Furthermore,  it  is  clear  that  □  is  an  inclusion  monotone 
interval  transformation,  and  that  the  interval  defined  by  (3.6)  is  the  element  of  IG 
of  minimum  width  which  contains  [a,bj.  If  x  £  RG  i*  a  real  number,  then 

(3.7)  Qc  =  [7x,AxJ 

will  be  called  the  HZpHeientation  of  x  in  IG  of  minimal  width.  If  [7x,Ax]  »  [c,d]  is 
nondegenerate,  then  [c,d]  also  represents  each  real  number  z  such  that  c  <  z  <  d  in 
this  way. 

Thus,  if  $  is  an  interval  transformation  which  maps  IRG  into  IRG,  then  the  cor¬ 
responding  noundid  tHont^oHiration 


w  *'■?**£ 


(3.8) 


maps  IRG  (and  hence  IG)  into  IG.  Interval  transformations  F:  D  C  IG  -*  IG  are  said  to 
be  computable  on  G;  furthermore,  if  F  has  property  (2.3)  (inclusion  of  the  real  trans¬ 
formation  $)  and  (2.5)  (inclusion  monotonicity),  then  F  will  be  called  a  computable 
■interval  extension  of  $ .  It  follows  that  one  can  start  with  an  interval  transformation 
♦  s  D  C  irg  -*  IRG  with  properties  (2.3)  and  (2.5)  (in  particular,  *  can  be  an  inclusion 
monotone  interval  extension  of  t) ,  and  obtain  the  computable  interval  extension  (3.8) 
of  p  by  applying  the  rounding  operator  D.  This  construction  can  be  carried  out  auto¬ 
matically  if  interval  arithmetic  and  a  library  of  computable  interval  extensions  of 
the  real  functions  entering  into  the  definition  of  i  are  available.  One  problem  which 
arises  in  this  connection  is  that  many  computers  lack  the  roundings  (3.4)  and  (3.5) 
required  to  support  interval  arithmetic.  This  defect  can  be  remedied  by  software 
[25],  [26]  or  microprogramming  [13],  if  necessary.  The  scientific  computing  language 
PASCAL-SC  [8]  includes  properly  rounded  interval  arithmetic  as  a  standard  feature. 
Henceforth,  all  interval  transformations  considered  will  be  assumed  to  be  computable. 

4.  Interval  iteration.  A  standard  approach  to  the  solution  of  equation  (1.1)  by  in¬ 
terval  methods  is  to  transform  it  into  a  fixed  point  problem  (1.2)  in  some  way,  and 
then  apply  inteAoaJt  iteration  [1],  [12],  This  process  takes  the  following  form: 

Given  an  interval  transformation  t  which  includes  i  in  the  sense  of  (2.3),  end  an  in¬ 
terval  XQ  thought  to  contain  a  fixed  point  x*  of  p,  the  sequence  of  intervals  {X^ }  de¬ 
fined  by 

(4.1)  X„.,  -  ),  n  -  0,1,2,..., 

n+i  n  n 

is  generated.  The  following  results  are  well  known  [24]: 

(i)  If  x*  e  XQ,  then 

00 

(4.2)  x*  e  x  -  n  xn, 

n-0 

so  that  existence  of  a  fixed  point  of  p  in  XQ  implies  convergence  of  (4.1)  in  the  sense 
that  X  ?  0,  the  empty  set. 

(ii)  If,  for  some  positive  integer  N, 

(4.3)  XN  -  0, 


TV  v  I 
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then  the  interval  XQ  contains  no  fixed  point  x*  of  4.  Hence,  divergence  of  (4.1)  in 
the  sense  of  (4.3)  implies  nonexistence  of  a  fixed  point  of  4  in  XQ;  this  is  the  con¬ 
trapositive  of  the  assertion  in  (i) . 

In  finite  dimensions,  one  has  further: 

(iii)  If 

(4.4)  4  (X  )  C  X„ 

for  some  finite  integer  M,  then  there  exists  a  fixed  point  x*  of  4  in  X  ,  and  (4.2) 

M 

holds  as  a  consequence  [10].  This  is  because  (4,4)  implies  that  the  continuous  opera¬ 
tor  0  maps  the  closed  compact  region  x  into  itself,  and  thus  x*  6  X  C  x_  by  the 

M  MO 

Schauder  fixed  point  theorem. 

Moreover,  if  x*  €  X  is  a  fixed  point  of  4,  and  x  «  [x  ,x  ] ,  X  “  [x,x] ,  then 
0  n  —n  n  — 

one  has  the  lower  and  upper  beinds , 

(4.5)  x  S  x*  S  x  ,  n  -  0,1,2,...,  and  x  £  x*  S  x, 

— n  n  — 

where  the  inequalities  are  interpreted  componentwise  in  several  dimensions.  From 

(4.5)  ,  approximations  to  x*  and  guaranteed  error  bounds  can  be  obtained  easily,  if 
desired. 

The  above  theory  requires  no  modification  in  case  the  calculations  are  performed 
on  a  finite  set  of  points  G  in  one  or  more  dimensions  by  taking  a  computable  interval 
extension  4  of  4 ,  and  the  initial  interval  X^  €  IG.  In  this  case,  however,  the  in¬ 
terval  iteration  process  (4.1)  is  iiruXi,  and  will  ZeAminate.  in  either  the  condition 
(4.3),  establishing  nonexistence  of  x*  in  X^,  or 


(4.6) 


X  =  x  a  0 
N+l  N 


in  which  case  X  =  X„  in  (4.2).  Condition  (4.6)  is  called  the  NickeX  ZoAminaXion  vU- 
N 

tQJii.on  [15].  The  upper  bound 


(4.7) 


N  <  G# (XQ) 


has  been  given  for  N  [24],  where  G#(XQ)  denotes  the  number  of  points  of  G  contained  in 


Interval  iteration  is  thus  an  obitAvabli  process  in  actual  computation,  since 
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it  yields  a  sequence  (xnln,0 


the  usefulness  of  which  can  often  be  determined 


by  inspection.  One  gets  nonexistence  of  x*  if  (4.1)  holds,  or  existence  if  (4.4)  is 
true  for  MSN,  and  the  guaranteed  bounds  (4.5)  if  XQ  contains  one  or  more  fixed  points 
x*  of  4.  One  is  in  the  dark  concerning  existence,  however,  if  (4.6)  holds  without 
(4.4)  occuring  along  the  way.  In  particular,  if  X^  C  ttX^)  properly,  then  X^  a  X^  =  X, 
and  the  interval  iteration  is  said  to  have  itaZttd.  The  difficulty  concerning  exis¬ 
tence  can  be  resolved  by  applying  seme  other  existence  test  than  (4.4)  to  X  =  X^,  a 
region  which  will  be  smaller  than  XQ  unless  the  iteration  stalled,  or  X  can  be  divided 
into  subregions  for  further  investigation,  using  the  algorithm  of  Moore  and  Jones, 
which  has  also  been  shown  to  be  finite  if  G  is  sufficiently  fine  [13),  [14). 

5.  An  interval  iteration  operator.  A  simple  way  to  convert  the  equation  (1.1) 
to  be  solved  into  a  fixed  point  problem  (1.2)  is  by  introduction  of  the  iteration  op- 
eACLton.  4  defined  by 

(5.1)  4  (x)  -  x  -  Yf  (x)  , 

where  Y  is  a  nonsingular  real  matrix,  in  which  case  the  two  problems  are  equivalent. 
Supposing  that  f  has  the  (Fr^chet)  derivative  f’,  represented  by  the  Jacobian  maVUx 
f '  (x)  =  Of^(x)/3Xj)  [20],  then  4  will  also  be  differentiable,  with 

(5.2)  4'  -  I  -  Yf • , 

where  I  denotes  the  identity  matrix.  The  derivative  (5.2)  can  be  used  to  construct 
the  mean  value  form  (2.7)  of  4,  provided  that  f'  has  an  inclusion  monotone  interval 
extension  F',  so  that  K  defined  by 

(5.3)  K(X)  -  y  -  Yf  (y)  +  [I  -  YF'(X)}*(X  -  y)  ,  y  »  m(X)  , 

will  be  an  inclusion  monotone  interval  extension  of  the  iteration  operator  4  in  (5.1). 

K  is  often  called  the  Kaawczifk  npeAatoA  [  5  ]  . 

An  alternative  form  of  the  Krawczyk  operator  (5.3)  can  be  obtained  by  the  use  of 
(5.1)  and  (2.15) .  One  has 

(5.4)  K (X)  *  4  (m (X) )  +  (I  -  YF*  (X)  )*w(X)  — j— 

which  expresses  K(X)  as  the  sum  of  the  ttZAatton  point  $(y)  of  y  *  m(X)  and  a  symmetric 
interval  vector.  Some  specific  choices  of  the  real  matrix  Y  will  now  be  examined. 
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(i)  An  ideal  choice  of  Y  is  Y  =  tf'(x*)3 


the  inverse  of  the  Jacobian  matrix 


-1 

of  f  at  a  solution  x  =  x*  of  equation  (1.1).  For  this  fixed  matrix,  the  ordinary  itera¬ 
tion  process  (1.3)  with  4  given  by  (5.1)  will  converge  quadratically  for  xQ  in  seme 
neighborhood  N(x*)  of  x*  [19).  This  gives  the  possiblity  of  quadratic  convergence  of 
the  endpoints  of  Xn  to  x»,  at  least  for  XQ  C  N(x*>.  Of  course,  since  x*  is  unknown,  it 
is  usually  not  possible  to  find  f'(x*).  However,  if  f’  satisfies  an  equation  of  the 
form  f'(x)  =■  h(f  (x) )  ,  then  one  has  f  •  (x*)  -  h(0),  and  V  can  be  found  without  knowing 
x*,  thus  giving  a  rapidly  convergent  iteration  method  of  simple  form  (191.  A  numerical 
example  of  this  type,  in  which  the  interval  iteration  converges  quadratically,  will  be 
given  below. 

(ii)  By  analogy  with  Newton's  method  [5  ),  one  can  choose,  if  it  exists, 

(5.5)  Y  -  If'  (y)f1  =  If'  (m(X))]'1. 

In  this  case,  the  iteration  point  <t>  (m (X> )  will  be  the  same  as  obtained  by  applying 
one  step  of  the  Newton  iteration  process  to  (1.1),  starting  from  xQ  ■  y  ■=  m(x)  .  The 
corresponding  interval  iteration  operator  K  defined  by 

(5.6)  K (X)  -  <Mm(X))  +  {I  -  [f'(B(X))j"1F'(X))'w(X)-i::i^- 

will  be  called  the  Ulutton  $OAjn  of  the  Krawczyk  iteration  operator  defined  by  (5.4). 

It  should  be  noted  that  the  discussion  of  interval  iteration  in  §4  also  applies  to 
the  case  in  which  Y  depends  on  X  as  long  as  0 (X)  contains  0(X)  for  an  operator  0  with 
fixed  point  x* ,  in  particular,  for  0  as  given  by  (5.1). 

(iii)  Another  possible  choice  of  Y  is 

(5.7)  Y  -  (rafFMX))]-1, 

again  provided  that  the  indicated  inverse  exists.  Writing 

(5.8)  F' (X)  »  m(F’ (X) )  ♦  jw(F' (X) ) •  1-1,1) 

by  (2.11),  substitution  into  (5.4)  yields,  with  .he  aid  of  (2.12), 

(5.9)  K  (X)  -  4  (m  (X)  )  +  |[m(F'(X)))  1  iw(F'  (X)  )w(X)  •  , 

where  |v|  *  (ly^l)  for  the  real  matrix  V  -  (y^)  .  The  interval  operator  K  defined 
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by  (5.9)  will  b«  called  the  Moow-Jonti  io>u n  of  the  Krawczyk  iteration  operator  [14]. 

Compared  to  (5.6),  (5.9)  has  the  computational  advantage  that  the  elements  of  m(P'(X)) 

can  be  obtained  by  simple  arithmetic  once  F‘ (X)  has  been  formed,  while  construction  of 

2 

f ' (m (X) )  requires  the  computation  of  the  v  elements  of  the  v*v  Jacobian  matrix  f'(y). 

(A  case  in  which  the  two  operators  coincide  will  be  considered  later.) 

The  Krawczyk  operator  is  useful  in  two  ways  in  connection  with  ordinary  iteration: 

(i)  K  may  be  used  to  ttit  an  interval  X  for  the  existence  or  nonexistence  of  a 
solution  x*  €  x  or  (1.1),  according  as  K(X)  C  x  or  K(X)  nX  •  0  [14].  In  case  this 

test  shows  that  x*  €  x,  this  establishes  the  suitability  of  xQ  ■  y  =  m(X)  as  a  start-  , 

ing  point  for  the  ordinary  iteration  (1.3)  under  fairly  mild  conditions  [14],  On  the 

I 

J  other  hand,  given  y  ■  x^,  one  may  take  X  =  X^  to  be  the  Cube,  with  diameter  2p  (in  the 

1*1  norm  for  R  )  centered  on  y, 

(5.10)  X  *  y  +  pe« (-1,1] , 

P 

for  which  (5.6)  and  (5.9)  become 

(5.11)  K(x  )  -  i(i  (y)  +  o{I  -  [f '  (y)  l”1*"  (X  )  [-1,1] « 

p  p 

where  e  *  (1,1,...,1),  and 

(5.12)  K(X  )  -  (My)  +  [m(F'  (X  ))]'1|W(F'  (X  >)e-  t-l,U  . 

p  <  p  p 

l 

respectively.  These  are  very  simple  to  compute,  since  the  formation  of  the  product 
Ae  for  a  real  or  interval  matrix  A  does  not  require  any  multiplications;  Ae  is  simply 
*  the  vector  consisting  of  the  sums  of  the  rows  of  A. 

If  one  wishes  to  examine  K(X  )  for  several  values  of  p,  then  the  Newton  form 

P 

(5.11)  has  the  computational  advantage  that  (f’(y)l  1  remains  fixed,  while  m(F'(Xp))] 

must  be  recalculated  in  the  Moore-Jones  form  (5.12)  for  each  p.  i 

(li)  K  may  be  used  to  determine  the  accuA&cy  of  the  elements  of  an  ordinary  | 

iteration  sequence  {x  }.  In  this  case,  for  y  *  ,  $ (y)  *  x  . .  A  choice  of  p  which  ■ 

n  n  n-*- a  •  j 

is  supported  by  the  Kantorovich  theory  for  the  convergence  of  Newton's  method  (21]  is  I 

M  -  2r  ,  where  j 

n  •  | 

<5*13)  *  1  xn  “  ♦IV1-  *  ,xn  "  *n+l,a>'  n  *  °'1'2 .  j 
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If  K (X  )  C  x  for  this  (or  other)  choice  of  p,  then  the  existence  of  x*  is  verified, 
and  one  has  the  error  bounds 

(5.14)  8x*  -  x  U  <  p,  (lx*  -  x  .1  S  -V (K (X  )), 

n  00  n+1  ®  2  p 

since  x*  €  K (X  )  C  x  .  From  (5.12), 

P 

(5.15)  "x*  -  xntl'»  S  §^|  tm(P*(Xp))l'1|w(F*(Xp))loi>, 

and  (5.11),  (5.12),  or,  more  generally,  (5.6),  also  give  componentwise  bounds  for 
x*  -  x 

n-»  1 

Thus,  the  interval  iteration  operator  K  can  be  ptumppLizd  to  scout  for  a  likely 
region  m  which  ordinary  iteration  will  converge  to  a  fixed  point,  or  po&ta.ppLied  to 
verify  the  existence  of  x*  and  obtain  error  bounds  for  a  result  obtained  by  ordinary 
iteration,  perhaps  computed  only  approximately.  In  this  respect,  the  Krawczyk  operator 
(5.3)  requires  a  slight  modification  for  computation  on  a  finite  set  G.  Suppose  that 
F  is  a  computable  interval  extension  of  f,  and  F'  is  a  computable  interval  extension 
of  f'.  Then,  one  can  use  K  defined  by 

(5.16)  K(X)  *  y  -  YF(y)  +  (I  -  YF’(X)}*(X  -  y) , 

where  all  arithmetic  operations  indicated  are  performed  in  correctly  rounded  interval 
arithmetic.  One  can  write  a  ffl  b  -  Q(a  ♦  b)  for  interval  addition  followed  by  directed 
rounding,  etc.  The  directed  rounding  operator  0  will  be  suppressed,  however,  for  clar¬ 
ity  of  notation,  and  its  use  will  be  understood  in  connection  with  actual  computation. 
The  computable  transformation  (5.16)  will  include  the  interval  transformation  (5.3), 
which  is  all  that  is  required  for  the  purpose  of  interval  iteration.  It  is  also  pos¬ 
sible  to  construct  a  computable  interval  extension  of  K  which  is  inclusion  monotone 
on  subsets  of  X^.  Suppose  that  for  each  y  €  xq,  one  can  find  z(y)  €  g  such  that 

(5.17)  #(y)  €  z (y)  +  ee*(-l,l] 
for  some  fixed  r  *  o.  Then,  from  (5.4), 

(5. IB)  K ( X)  *  z (y)  +  ({I  -  YF' (X) )*|w(X)  +  te)*[-l.l] 

is  a  computable  interval  extension  of  the  Krawczyk  operator  which  has  the  property 
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x  c  z  c  xQ  “  K (X)  C  k(Z)  ,  which  will  be  called  inclusion  monotorUciXy  on  iubieti  oi 
XQ.  Since  interval  iteration  stays  entirely  within  the  initial  interval  XQ,  this  is 
satisfactory  for  practical  purposes. 

For  the  Moore-Jones  transformation  (5.12),  this  enlargement  of  the  Krawczyk  trans¬ 
formation  takes  the  simple  form 

(5.19)  K(xp)  =  z(y>  +  {el  +  ||y|w(F' (Xp)) ie- [-1,1) 

where  [m(F' (X  ))!  1  £  Y.  Thus,  once  F' (X  )  and  Y  have  been  obtained,  the  only  inter- 

P  P 

val  operation  to  be  performed  in  the  evaluation  of  (5.19)  is  the  trivial  multiplica¬ 
tion  by  (-1,1J;  the  rest  can  be  done  with  properly  rounded  real  arithmetic.  In  par¬ 
ticular,  calculation  of  the  components  of  the  vector  (el  +  y|w(f' (X^) ) }e  requires 
only  operations  with  non-negative  numbers  and  upward  rounding  (3.4).  Recall  that  |y| 
is  defined  componentwise  by  (2.9). 

6.  Interval  iteration  with  fixed  Y.  What  will  be  considered  here  is  interval 
iteration  with  the  simple  case  (i)  of  the  Krawczyk  transformation  defined  in  55,  in 
which  the  matrix  Y  remains  fixed.  One  would  like,  of  course,  to  have  Y  »  lf'(x*)]  \ 
or  an  approximation  to  the  inverse  of  the  Jacobian  matrix  of  f  at  x  *  x*,  <  solution 
of  equation  (1.1).  This  is  possible  in  the  special  case  that  (1.1)  is  a  Lineal  sys¬ 
tem  of  equations,  that  is,  f  is  an  a^ine  operator  given  by 

(6.1)  ftx)  -  Ax  -  c, 

where  A  is  a  nonsingular  matrix  and  c  is  a  given  vector.  Here, 

(6.2)  f  (x)  -  A, 

a  constant  operator  [20],  and  thus  [f'(x*)l  1  «  A  *.  With  this  choice  of  Y,  it  follows 
from  (5.3)  that  K(Xg)  »  x»,  and  thus  the  interval  iteration  will  converge  to  the  solu¬ 
tion  x  -  x*  of  Ax  -  c  in  one  step  for  arbitrary  X^  such  that  x*  £  X^.  In  practice, 

Y  will  only  be  an  approximation  to  A  Because  of  (6.2),  one  can  take  F'(X)  -  A  for 
arbitrary  X,  and  thus  (5.4)  becomes 

(6.3)  K (X)  -  m(X)  -  Y  (Am(X)  -  c)  +  {I  -  YA}w(X)  •— 

or,  in  case  X  is  taken  to  be  a  cube  X  -  y  +  pe*(-l,l),  one  has 

o 
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(6.4) 


K(X  )  =  y  -  Y (Ay  -  c)  +  p j I  -  YA|e-[-l,l], 
which  is  a  very  simple  interval  transformation  to  calculate.  The  XtiiAue.  r  =  f(y)  = 
Ay  -  c  enters  into  (6.4)  in  a  natural  way. 

In  connection  with  the  solution  of  linear  equations,  it  should  be  noted  that  in¬ 
tervals  do  not  form  a  linear  space.  From  the  rule  for  iubfiaction  [9),  Ill), 

(6.5)  [a,b]  -  [c,d]  =  [a  -  d,  b  -  o], 

it  follows,  for  example,  that 

(6.6)  [0,1]  -  [0,1]  -  [-1,1]. 

Thus,  techniques  from  algebra  which  depend  on  linearity  cannot  be  expected  to  extend 
to  interval  analysis,  since  a  linear  substructure  is  not  present,  sane  failures  in 
attempted  solution  of  linear  systems  by  interval  versions  of  methods  from  linear  alge 
bra  are  due  to  a  lack  of  understanding  of  this  fact,  not  to  a  defect  of  interval  anal 
ysis. 


For  an  example  of  a  simple  system  of  nonlinear  equations  to  which  the  method  of 
this  section  applies,  consider  the  equations 


(6.7) 


exp{-u  +  v)  -  0.1  ■  0, 
exp{-u  -  v)  -  0.1  “  0, 


due  to  Cuyt  and  Van  der  Cruyssen  [  4]  for  x  -  (u»v)  €  R2.  (Perhaps  this  system  is 
too  simple;  an  equivalent  linear  system  can  be  obtained  by  taking  logarithms.)  The 
Jacobian  matrix  corresponding  to  (6.7)  is 


(6.0) 

with  inverse 


J (u,v) 


C-  exp{-u  +  v) 
-  exp{-u  -  v) 


exp{-u  +  v) 
-  espl-u  -  v) 


) 


(6.9) 


J(u,v) 


exp{u  -  v] 
exp(u  -  v) 


-  exp{u  + 

-  exp{u  ♦ 


Thus,  if  x*  »  (u*,v»)  satisfies  (6.7),  then 
(6.10)  Y  -  J(u*,v*)_1  - 
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is  obtained  without  knowing  x*.  Using  (6.10),  the  iteration  process  (1.3)  for  <J>  given 
by  (5.1)  with  f  defined  by  (6.7)  is 


u  ,  -  u  ♦  5exp(-u  }(exp{v  }  +  exp{-v  })  -  1, 

(6.11)  n+l  n  "  n 

v  ,  ■  v  -  5exp(-u  )(exp{v  )  -  exp{-v  )), 
n+l  n  n  n  n 

which  will  converge  quadratically  if  xQ  ■  (Ug,Vg)  *s  c*OS€  ®nou9h  to  x#  *  (u*,v*)  (19]. 
The  c)ioico  xQ  *  (4.3  ,2.0)  t  4  )  leads  to  the  results  shown  in  Table  6.1.  The  calcula¬ 
tions  were  done  with  an  HP-33E  pocket  calculator,  with  rounding  to  four  decimal  places. 


r — 

n 

u 

n 

V 

n 

n  ■  lx  .  -  x  1 
n  n+l  n  ® 

0 

4 . 3000 

2.0000 

0.4921 

1 

3.8105 

1.5079 

0.4755 

2 

3.3350 

1.0324 

0.4366 

3 

2.8984 

0.5958 

0.3481 

4 

2.5503 

0.2477 

0.1954 

5 

2.3549 

0.0524 

0.0498 

6 

2.3052 

0.0026 

0.0026 

7 

2.3026 

6.9985  (-06) 

- 

8 

2 . 3026 

-4.6502  (-10) 

- 

9 

2 . 3026 

2.8499  (-10) 

- 

10 

2 . 3026 

2.8499  (-10) 

- 

Table  6.1.  Numerical  Results  for  the  Iteration  (6.11). 

The  elements  of  the  iteration  sequence  (x  }  can  be  used  to  form  cubes 

n 

(6.12)  X  -  (U,V)  -  x  +  2n  e*  [-1,1]  -  (u  ±  2n  ,  v  ±  2n  )  , 

n  n  n  n  n  n 

which  have  the  Krawczyk  transformations 

(6.13)  K (X)  -  x  .  +  n  (I  -  YJ(U,V) }e* (-1,1) , 

n+i  n 

using  (5.4).  Interval  iteration  starting  with  (6.12)  as  XQ  and  using  (6.13)  will  stall 


for  n  =  0,1 ,2,3,4.  For  n  =  5,  however,  one  gets 

16.14)  XQ  D  K(Xq)  *  vi  3  K(Xx)  *  X2  D  K(X2)  =  X3, 

with  '<  *  (U  ,V  )  given  in  Table  6.2. 

n  n  n 


Table  6.2.  Results  of  Interval  Iteration. 


The  results  given  here  verify  the  existence  of  x*  e  x  ,  n  ■  0,1, 2, 3,  and  the 

approximate  solution  x  -  (2.3026  ,0.0000)  is  guaranteed  to  be  accurate  to  four  deci¬ 

mal  places,  since  x  *  m(X3)  ,  jw(X3>  *  0.00002.  The  calculations  here  were  also  done 
using  an  HP-33E,  with  the  subroutines  [23)  augmented  by  simple  programs  for  exp{X}  and 
exp{-x},  and,  of  course,  directed  rounding  of  the  results  to  the  number  of  places  in- 
uicated  in  the  table. 

7.  The  eigenvalue-eigenvector  problem.  Interval  iteration  can  also  be  applied 
to  the  eigenvalue -zigenvectoX  pxobl&n  for  a  real  matrix  A,  which  is  to  find  real  num¬ 
bers  *  and  corresponding  nonzero  vectors  x  such  that 
(7.1)  Ax  -  >x  -  0. 

The  condition  x  ¥  o  is  often  enforced  by  a  noX/naCizdCion  of  x, 

(7.^)  (a)  (x,x)  -  1  -  0,  or  (b)  (a,x)  -  1  «  0, 

the  Latter  for  a  fixed  vector  a.  The  system  of  equations  formed  by  (7.1)  and  (7.2) (a) 

or  (7.2)  (b)  is  a  nonlinear  system  of  v  +  1  equations  in  the  v  ♦  1  unknowns  x^ ,  , . . . , 

x  ,  * .  In  fact,  the  result  is  a  special  case  of  what  are  called  quadratic  systems  of 
equations,  in  which  the  transformation  f(x)  has  the  form 
( 7 . j)  f (x)  »  Bxx  +  Lx  +  c, 

in  which  B  «  (b_k)  is  a  bitinQOA  vpeAatoX ,  L  is  a  given  matrix,  and  c  is  a  fixed 
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vector  [17]  ,  [20].  The  dzfu.va.ti.vi  of  the  quadfiatiz  opihJJltoK  f  defined  by  (7.3)  is 

(7.4)  f • (x)  -  2Bx  +  L, 

where  B  =*  j(b_k  +  bikj)  is  a  4 yrmitsUc  bilinear  operator  [17],  [20].  It  follows  that 

(7.5)  F* (X)  -  2B-X  +  L 

gives  an  interval  extension  of  (7.4),  furthermore,  since  B  is  a  linear  operator  on  x 
(into  the  space  of  matrices) ,  one  has 

(7.6)  m(F'(X))  ■  2Bm(X)  +  L  ■  F'  (m(X) }  , 

which  means  that  the  Newton  form  (5.6)  and  the  Moore-Jones  form  (5.9)  of  the  Krawczyk 
iteration  operator  will  coincide  for  quadratic  systems  of  equations.  Thus,  an  inter¬ 
val  iteration  method  based  on  Newton's  method  will  have  a  simple  computational  form  in 
this  case.  In  particular,  the  formulation  of  the  eigenvalue-eigenvector  problem  (7.1) 
as  a  system  of  quadratic  equations  [2]  ,  [18]  leads  to  a  simple  interval  iteration 
method,  which  will  be  given  below. 

Returning  to  (7.1),  if  x*  ?  0  is  an  eigenvector  of  A,  then  the  corresponding  eigen¬ 
value  X*  is  given  by  the  R ayteXgh  quotiint 

(7.7)  X*  -  R(x*)  -  (Ax*,x*)/(x*,x») . 

Consequently,  if  x*  6  X  for  some  interval  vector  X  such  that  0  6  X,  then  the  inttfivaZ 
RayleXgh  quotient 

(7.8)  R (X)  -  (AX,X)/(X,X) 

can  be  used  to  obtain  louHK  and  upptA  boundi  for  X*:  Since  X*  6  R(X)  •  [c,d] ,  it  fol¬ 
lows  that 

(7.9)  c  s  X*  s  d. 

R (X)  can  be  computed  automatically,  using  interval  arithmetic,  with  the  extension 

r  2 

(7.10)  (X,X)  *  l  X 

i-1  1 

of  (x,x)  ,  where  the  iquaAl  of  an  interval  is  defined  to  be 

,  [min{a2 ,b2} ,max{a2 ,b2) ]  if  ab  >  0, 

(7.11)  [a,b]  = 

[O.maxta  ,b  )]  if  ab  <  0. 


-  16 


2 

Formula  (7.11)  gives  a  better  extension  of  f(x)  ■  x  than  F(X)  •  X»X  for  intervals  X 
which  contain  both  positive  and  negative  numbers;  for  example,  [-1,1]^  -  [0,1],  while 
[-1,11 •  [-1,11  -  1-1,1). 

The  eigenvalue-eigenvector  problem  for  A  will  be  considered  in  the  form  of  the 
quadratic  system  (7.1) - (7.2) (b) .  Eigenvectors  orthogonal  to  a  will  not  satisfy  this 
systen,  hut  can  be  obtained  from  a  system  with  a  different  choice  of  a.  It  will  be 
shown  that  the  Krawczyk  transformation  required  to  apply  interval  iteration  to  this 
problem  can  be  expressed  as  the  sum  of  the  point  obtained  by  Newton's  method  [  2 1 ,  [18] , 
and  an  easily  computed  symmetric  interval  vector. 

Following  [igl,  a  will  be  taken  to  be  a  unit  vector  in  (7.21(b),  for  example, 
a  -  (0,...,0,1).  This  gives  x^  •  1,  and  (7.1)-(7.2) (b)  becomes  the  system  of  v  equa¬ 
tions 


(“ll  -  X)X1  +  #12X2 


..  + 


1, v-1  v-1 


+  a 


lv 


0, 


a21Xl  +  **22  *  X,X2  +  •••  ♦ 


2, v-1  v-1 


+  a. 


2v 


0, 


(7.12) 


a  ,  +  a  ,  ,x,  +  ...  +  (a  .  ,  -  A)x  ,  +  a  ,  -  0, 

v-1,1  1  v-1, 2  2  v-1, v-1  v-1  v-l,v 


+  ...  + 


a  ,x  . 
v,v-l  v-1 


+  a  -  A  »  0, 

vv 


for  the  v  unknowns  x^  x^ ,  ...,  x^_^ ,  A.  The  derivative  f’  of  the  operator  f  defined 
by  the  left  side  of  (7.12)  is 


(7.13) 


r 


a, ,  -  A 


t’(V 


12 

1 ,  v— 1 

i 

*21 

*22  *  X 

*2 ,v-l 

i 

X 

KJ 

“v— 1,1 

*v-l,2 

...  a  ,  .  -  A 

v-1, v-1 

v-1 

avl 

“v2 

...  a  , 

v,v- 1 

-1 

where  C  «  (n^,. . . ,x  , A) .  An  Interval  extension  F' (=)  of  (7.13)  is  obtained  by  simply 


replacing  r,  by  the  interval  vector  =  «  (x^,.... 


X^_^,A).  Since  most  of  the  components 


of  F*  (  =  )  are  real  lumbers,  the  real  matrix  w(F'(S))  has  the  simple  form 


w(A) 

0 

0 

w  (A) 

...  0 

...  0 

w(Xx) 

w(X2) 

0 

0 

...  W(A) 

W(X  .) 

V-1 

0 

0 

...  0 

1 

Thus,  the  vector  w(F‘(5))w(E)  is  given  by 


(7.15)  w(F'(S))w(E)  -  2w(A)  (w(X,i .  '  j-11  ,0)  “  2w  (A)w  (X (v) )  , 

where  X*v^  »  (X^ , . . . ,Xy_^ ,0) •  For  y  «  [f*  •>  n  “  m(E) ,  the  transformation 

(5.9)  becomes 

(7.16)  K(H)  -  n  -  yf  (n)  +  ;v(X<V*)  •  (-1.1J , 


the  sum  of  the  Neioton  pO-int  t  ■  n  *  *f(n)  and  a  simply  computed  symmetric  interval 
vector.  Formula  (7.16)  simplifies  further  for  wfXj)  «  w(Xj)  ■>...»  w(Xv-1>  “  2P< 
one  has 


(7.17)  K(E)  -  n  -  xf  (n)  +  pw(A)  | y  | e <v) » J-1,1) , 

where  e*'J'  «  (1,...,1,0). 

An  alternative  method  for  obtaining  lower  and  upper  bounds  for  eigenvalues  of  A 
is  available  from  (7.12).  If  x*  6  x  ■  (X^ , . . . ,X^_^ ,1)  is  an  eigenvector  of  A,  then 
A*  €  A  for  the  corresponding  eigenvalue  A*,  where 


v-1 

<7.ie>  A  -  a  +  T  a  •X. 

vv  ^  vj  j 

is  the  interval  version  of  the  last  equation  of  (7.12).  The  use  of  this  expression 
and  the  interval  Rayleigh  quotient  R(x)  given  by  (7.8)  to  obtain  bounds  for  eigen¬ 
values  of  matrices  will  be  illustrated  in  the  example  below. 

Consider  the  simple  example 


(7.19) 


:  j! 
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Here#  system  (7.12)  ,  written  for  £  =  (u,v,A),  becomes 
3u  +  v-.*u-l  =  Q, 

(7.20)  u  +  5v  -  Av  -  1  =*  0, 

-u  -  v  “  A  +  3-0. 


For  the  initial  approximations  uQ  *  vQ  =  1,  the  third  equation  of  (7.20)  gives  A  ■  1. 
The  results  of  applying  Newton's  method  to  (7.20),  starting  from  *  (1,1,1),  are 
shown  in  Table  7.1.  These  are  rounded  to  four  places  from  a  double -precision  calcu¬ 
lation  performed  on  a  UNIVAC  1160,  using  the  program  NEWTON  [  6). 


n 

u 

n 

V 

n 

A 

n 

0 

1.0000 

1.0000 

1.0000 

i 

0.8182 

0.2727 

1.9091 

2 

0.9755 

0.0166 

2.0079 

3 

1.0002 

-0.0001 

1.9999 

4 

1.0000 

0.0000 

2.0000 

Table  7.1.  Results  of  Newton's  Method. 

In  order  to  construct  an  interval  E  to  which  to  apply  the  Krawczyk  transformation 
(7.17),  it  is  helpful  to  split  the  vector  £  into  its  vtct0>l-£'ikz  part  ■  (*1,X2' 

••’,Xv-l'0),  3X1(3  *ts  va£ue-£4fee  part  Cv  *  (0,...r0,A)  -  Ae^,  *  (0,...,0,1).  Then, 
for 


(7.21) 


n  *=  |r(v>  -  f(v)l 
n0  ~  ^n+1  n 


it  is  reasonable  to  choose  p  *  2^#  w(A)  ■  4|a^^  ”  *nl*  This  9i-ves 


(7.22) 


:  -  ?n  +  2nQe(v) -1-1,1]  +  t-1 ,1) 


as  the  interval  to  be  examined.  Splitting  the  vector 

(7.23)  t  =  I  y  I  e  <v>  -  I  [m(F' (H)]'1|e<V>  -  t(v)  *  -  t<V>  *  Te^, 

it  follows  from  (7.17)  that  the  condition  k(5)  C  5  will  hold  if  the  inequalities 


n051/2T'  1\,+1  *  V  4  1/8,t 
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(7.24) 


are  satisfied.  Observation  of  the  sequence  given  in  Table  7.1  shows  that  (7.24)  holds 


for  n  =  1 ,  thus ,  for 


(7.25) 


one  has 


(7.26) 


.(3) 


A[0 . 3060,1. 3304] \ 
[-.2395,0.7049]  j, 
V.C1. 5139, 2. 3043]/ 

[0.7469,1.2041]\ 
[-.1671,0.2003]  ]  C  E. 
[1.0710,2.1448]/ 


The  non-negative  numbers  entering  into  the  calculation  have,  of  course,  been  rounded 
upward  to  four  decimal  places,  so  that  (7.26)  proves  that  the  matrix  A  given  by  (7.19) 
has  an  eigenvector  x*  «  (xj ,x* ,1)  and  a  corresponding  eigenvalue  X*  which  satisfy 

(7.27)  0.7469  <  x*  S  1.2041,  -.1671  S  X*  S  0.2003,  1.8710  S  X*  S  2.1448, 
since  £;*  €  K ( H)  ,  and  the  interval  iteration  using  the  transformation  (7.17)  will  con¬ 
verge.  The  interval  Rayleigh  quotient  applied  to  X  -  ( [0.7469,1 .2041] ,  [-.1671 ,0.2003] , 
1)  obtained  from  K(E)  in  (7.26)  gives 

(7.28)  R(X)  -  [0.8752,3.8826], 


while  formula  (7.18)  gives 

(7.29)  ^  -  [1.5956,2.4202] . 

The  bounds  that  (7.28)  and  (7.29)  provide  for  X*  are  rigorous,  but  inferior  to  (7.27) 
in  this  case. 

The  transformation  (7.17)  can  also  be  used  to  check  the  final  result  L.  of  the 

4 

numerical  computation  given  in  Table  7.1.  Take,  for  example, 

(7.30)  5  »  +  O.OOOle* [-1,1] , 

so  that  o  *  0.0001,  w(A)  =  0.0002.  Since  |y|e*3*  »  e  for  these  values,  one  has 
<7-31>  K(5>  -  +  0.0000002e- [-1 ,1]  C  =, 

which  proves  that  the  components  of  0  are  accurate  to  six  decimal  places.  Of  course, 
the  numerical  iteration  process  stumbled  onto  the  exact  solution  ■  £  in  this  case, 
but  this  fact  is  not  required  for  the  rigorous  error  analysis  (7.31). 
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